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Magnetic phase diagram of the S=l/2 antiferromagnetic zigzag spin chain in the 
strongly frustrated region: cusp and plateau 
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We determine the magnetic phase diagram of the S = 1/2 antiferromagnetic zigzag spin 
chain in the strongly frustrated region, using the density matrix renormalization group method. 
We find the magnetization plateau at 1/3 of the full moment accompanying the spontaneous 
symmetry breaking of the translation, the cusp singularities above and/or below the plateau, 
and the even-odd effect in the magnetization curve. We also discuss the formation mechanisms 
of the plateau and cusps briefly. 
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For the purpose of clarifying the role of the frustration 
in low-dimensional quantum spin systems, the S = 1/2 
antiferromagnetic zigzag spin chain has been attracting 
considerable attention, since it minimally contains the 
frustrating interaction without loss of the translational 
invariance. 1-4 ) The Hamiltonian of the model is given 

by 

H = 



[JxSi ■ S i+ i + hSi ■ S l+2 ] -H^2S?,(1) 



where S is the 5=1/2 spin operator, H is the mag- 
netic field, and J\ and J 2 denote the nearest and next 
nearest neighbor couplings, respectively. We introduce 
the notation a — J2/J1 for simplicity. Recently, the 
zigzag chain was realized as SrCuC^, 5 -' Cu(ampy)Br2, 6 * ) 
(N 2 H 5 )CuCl 3 7 ) and F 2 PIMNH. 8 ) 

The Hamiltonian (1) has a very simple form, but cap- 
tures a variety of behaviors induced by the frustration. 
In particular, the zigzag chain in a magnetic field has 
been studied actively 9-13 ^ and it has been clarified that 
cusp singularities appear near the saturation field and/or 
in the low field region for a < 0.6, in accordance with the 
frustration-driven shape change of the dispersion curve of 
the elementary excitations. 12 ~ 14 ) However, the magnetic 
phase diagram including the strongly frustrated region 
(a > 0.6) has not been acquired yet. Here we remark 
that the magnetization curve of a = 0.6 is still quite 
different from the one in the a — > 00 limit (see Fig. 2 
(a)). Thus, as a is increased beyond a = 0.6, the in- 
trinsic structural change of the magnetization curve can 
be expected particularly around a ~ 1, where the most 
significant competition is achieved. 

In this paper we address the magnetization process of 
the zigzag chain in the strongly frustrated region(a > 
0.6), using the density matrix renormalization group 
(DMRG) method. 15) We then find that the obtained 
phase diagram exhibits rich physics, as is shown in Fig.l. 
For 0.56 < a < 1.25, the magnetization plateau appears 
at 1/3 of the full moment, accompanying the sponta- 
neous breaking of the translational symmetry with the 
period three. Moreover, the cusp singularities in the 
magnetization curve show quite interesting behavior; as 



a is increased, the high field cusp merges into the 1/3 
plateau at a ~ 0.82. Also the low field cusp merges 
the 1/3 plateau at a ~ 0.7, but it appears again when 
a > 0.7. In addition, we find an interesting even-odd 
effect in the magnetization curve for a > 0.7. 
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Fig. 1. Magnetic phase diagram of the zigzag spin chain. D: dimer 
gapped phase, P:l/3 plateau, TL1: one component TL liquid, 
TL2: two component TL liquid, and EO: even-odd behavior 
branch. The open circles denote the position of the cusp singu- 
larities. The solid circles represent the boundary of the dimer 
gapped phase. The triangles represent the upper and lower edges 
of the 1/3 plateau, and the crosses indicate the upper edge of 
the even-odd effect in the magnetization curve. The saturation 
field is normalized to be unity. 

In analyzing the phase diagram, an important feature 
of the zigzag chain is that the Hamiltonian (1) interpo- 
lates between the single Heisenberg chain(J2 = 0) and 
the double Heisenberg one(Ji = 0) continuously. In 
other words, we can describe the frustration effect in 
the zigzag chain as the interplay between the single and 
double chain natures. Although such a view point has 
not been emphasized so far, it provides some essential 
insights on the frustration effect in the magnetization 
curve. In fact, the cusp transition for a < 0.7 mentioned 
above is described as the transition between the one and 
two component Tomonaga-Luttinger(TL) liquid; 12 -* the 
two component TL liquid is realized below the lower cusp 
and/or above the upper cusp reflecting the two chain 
nature of the system, while the middle-field branch still 
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consists of the one component TL liquid. In the follow- 
ing, we show the results of the magnetization curves for 
a > 0.6 and then proceed to a detailed analysis of each 
phase. How can we connect the above noted competition 
with the characteristic behaviors in each phase? 

magnetization curves: In order to compute the magne- 
tization curve, we employ the DMRG of both the finite 
system size and infinite system size versions. However, 
here we mainly present the results for the finite size sys- 
tems. We compute the minimum energy En(M) for a 
magnetization M of the N = 192 site system with the 
free boundary condition, and then obtain the magnetiza- 
tion curve by determining the level crossing point from 
En{M) and En{M ±1). The maximum number of the 
retained bases in the DMRG computation is typically 60, 
where the energy is converged sufficiently. 

In Fig. 2 (a) we show the magnetization curve for 
a = 0.6, where the scale of the magnetic field is normal- 
ized with the saturation field H s . The cusp singularities 
at H/H s ~ 0.22 and H/H s ~ 0.75 are the same as those 
reported in the previous papers, which are explained well 
by the shape change of the spinon or spin wave disper- 
sion curves. 13 ' 14 ^ A novel point in the figure is that the 
1/3 plateau is found to appear accompanying the sponta- 
neous symmetry breaking of the translational invariance. 
The lower critical field of this plateau is H c \/H s = 0.424 
and the upper critical field is H C 2/H S = 0.477. The de- 
tails of the 1/3 plateau are discussed later. 

As a is increased, the cusps approach the plateau, 
which implies that the two component TL region be- 
comes more dominant, and at the same time the width 
of the plateau extends. Figure 2 (b) shows the magne- 
tization curve of a = 0.7, in which the width of the 1/3 
plateau is extended to 0.209. An anomalous step cor- 
responding to M — 33 in the 1/3 plateau is due to the 
open boundary effect. By analyzing the size dependence, 
we have confirmed that this step is almost independent 
of N and is thus negligible in the bulk limit. In addi- 
tion to the plateau, an interesting point is that the low 
field cusp is merged into the 1/3 plateau, and then the 
branch below the plateau sticks into the plateau with in- 
finite slope. On the other hand, the high field cusp still 
survives. 

In Fig. 2 (c), we show the result for a — 0.8. In the 
figure we can clearly see that the high field cusp almost 
merges into the 1/3 plateau; a precise computation indi- 
cates that the disappearing point of the high field cusp 
is a ~ 0.82. On the other hand, a cusp singularity ap- 
pears again in the low field region and, moreover, the 
magnetization curve above this low field cusp remark- 
ably exhibits an anomalous even-odd oscillation with re- 
spect to the magnetization M . Of course, the anomalous 
steps corresponding to M=odd must be skipped in the 
true magnetization curve of the finite size system which 
is determined from Ejq(M) with M =even. 16 ) However, 
we purposely show these anomalous odd steps, since they 
are reflecting an interesting aspect of the double-chains 
nature of the system. We discuss the details of this even- 
odd effect later, but here simply note that the magneti- 
zation curve obtained with the infinite system DMRG is 
consistent with the present 192-site result. 




Fig. 2. Magnetization curves for a = 0.6(a), 0.7(b) and 0.8(c). 
The vertical axis indicates m = M/N. 



As a is increased beyond a — 0.82, the 1/3 plateau 
turns to shrink. We show the results for a = 1.0 in 
Fig. 3 (a), where the plateau becomes narrower compared 
with that for the a = 0.8 result and the cusp shifts to 
the lower field side. However, the high field cusp does 
not appear again. Instead, the even-odd behavior of the 
magnetization curve is extended above the 1/3 plateau. 
This upper edge of the even-odd effect shifts to the higher 
field side as a is increased. Although there is a possibility 
of the cusp appearing at the upper edge of the even- 
odd branch, we can not confirm it within the present 
numerical results for 192 sites. 

Figure 3(b) shows the magnetization curve of a -1 = 
0.8. In the figure, the region of the even-odd oscilla- 
tion is extended to a higher field region and the plateau 
becomes indistinct. Indeed, for aT 1 = 0.8, the infinite 
system DMRG illustrates a continuous curve above the 
low filed cusp up to the saturation field, implying that 
there is no 1/3 plateau. As a is increased further, the 
branch above the cusp becomes more similar to that of 
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Fig. 3. Magnetization curves for a = 1.0(a) and a 1 = 0.8(b). 



the (double) Heisenberg chain. At the same time, the 
low field cusp itself approaches the dimer gap bound- 
ary. In the Ji — ► 0(a _1 — > 0) limit, the "amplitude" 
of the even-odd oscillation becomes smaller and finally 
the magnetization curve becomes equivalent to that of 
the Heisenberg chain again. However, we note that it is 
difficult to see how the cusp behaves in this limit simply 
by the direct observation of the magnetization curve. 

1/3 plateau state: Let us discuss the nature of the 1/3 
plateau. According to the quantization condition of the 
magnetization plateau, 17 ) the 1/3 plateau state satisfies 
(7(1/2— M) ^integer with M = 1/6, where q is the period 
of the plateau state. This implies that the translational 
symmetry of the Hamiltonian (1) must be spontaneously 
broken with q = 3, 6, ■ • • . In order to see it, we calculate 
the local spin moments at the plateau along the chain 
direction. The result is shown in Fig. 4, where the value 
of the local moment is developed to about (5?) ~ 0.38 
or —0.26 and thus the "up-up-down" type long-range 
order is realized. Moreover, the boundary effect decays 
rapidly, which supports that the 1/3 plateau state has 
the excitation gap. 

What is the origin of this up-up-down structure? We 
have also investigated the magnetization curve of the 
zigzag XXZ chain with the Ising-like anisotropy, and 
then found that the 1/3 plateau is extended up to the 
Ising limit. Here it should be remarked that the 1/3 
plateau of the up-up-down array is proven for the Ising 
anisotropic limit, 18 ) where the up- up-down structure is 
especially compatible with the triangle structure of the 
zigzag chain. In addition, we have also checked that a 
similar 1/3 plateau of the up-up-down array appears for 
the Ising-like classical XXZ zigzag chain reflecting the 



triangular structure. 19 - 1 However, this "classical" plateau 
vanishes at the isotropic XXX chain. Thus the present 
1/3 plateau is induced by the quantum effect. 




Fig. 4. The distribution of the local spin moment at M=32(l/3 
plateau) for a = 0.8 

even-odd effect: The anomalous even-odd behavior 
for a > 0.7 can be explained in terms of the forma- 
tion of the bound state of magnetic quasi particles as 
follows. The spin chain with a finite magnetization is 
described by the quasi particles carrying S z — 1 dis- 
tributed on the chain. 20 ) When a is small and the 
single chain nature is dominant, the magnetic particle 
is distributed on the single chain, but if a becomes 
large and the double chain nature becomes significant, 
the particles should be distributed on each chain equiv- 
alently. However, if the total magnetization is odd, 
the remaining one particle is sitting on one of the two 
chains. Thus the M=odd system is less stable than 
the M=even one, which suggests the formation of the 
bound state in the thermodynamic limit. In order to 
see it, we define the binding energy of the TV-site sys- 
tem as A N (m) ee 2(E N (M + 1) - E N (M)) - (E N (M + 
2) - E N (M)) = 2E N (M + 1) - E N (M) - E N (M + 2) 
for M =even. Figure 5 shows the obtained result with 
the size extrapolation of Atv(to) ~ A(m) + const /N 
for N = 240, 192, 168, 120, 96, 84, 21 ) which shows clearly 
that the binding energy develops in the region corre- 
sponding to the even-odd behavior. Here we note that 
this size extrapolation does not work well around the low 
field cusp(m ~ 0.6) due to the singular behavior. 

The formation of the bound state of the quasi parti- 
cles provides an essential insight into the cusp mecha- 
nism for a > 0.7. The magnetization curve below the 
low field cusp consists of the two-component TL liquid 
that is continuously connected from the a < 0.7 region, 
while the branch above the low field cusp consists of the 
bound state of the magnetic particles. Thus the cusp 
mechanism for a > 0.7 is described by the transition 
from the two component TL liquid to the single com- 
ponent TL liquid consisting of the bound state, which 
is completely different from the mechanism for a < 0.7 
based on the shape change of the dispersion curve of the 
spinon excitation carrying S z = 1/2. 

To summarize, we have presented the magnetic phase 
diagram of the 5=1/2 zigzag spin chain including the 
strongly frustrated region, based on the 192-site results. 
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Fig. 5. The binding energy for a = 1.0 obtained by the size extrap- 
olation: Ajv(m) ~ A(m)+const/A''. The open circles denote the 
result for A 7 =240 and 120. The upward and downward triangles 
denote those for N =192 and 96, and 168 and 84 respectively. 
These results are in good agreements with each other. 



Although we can not eliminate the possibility of a "fine 
structure" appearing in the magnetization curve beyond 
the resolution of the present calculation, we believe that 
the phase boundaries of the cusps and plateau are cor- 
rectly obtained. The resulting phase diagram in Fig.l 
exhibits a variety of interesting physics: The 1/3 plateau 
accompanying the spontaneous breaking of the transla- 
tional symmetry, the characteristic move of the cusp sin- 
gularities around the plateau and the even-odd behavior 
of the magnetization curve. In particular we have clar- 
ified two types of cusp mechanisms: the cusp based on 
the shape change of the dispersion curve and the cusp 
attributed to the formation of the bound state related 
to the even-odd behavior of the magnetization curve. In 
the background of these behaviors, we have focused on 
the role of the frustration in terms of the competition 
between the single and double chain properties. The 
magnetic field triggers the switching of these two na- 
tures which induces the cusp transitions. In addition, 
the 1/3 plateau and the cusps remarkably merge with 
each other at a ~ 0.7 and 0.82. However, the details 
of the mechanism for this remain an interesting future 
problem. 

Although the zigzag chain is a quite simple model, we 
consider that it captures the intrinsic aspects common to 



a class of frustrated spin chains having the zigzag type 
structure. For the exploration of such models, we believe 
that our result is instructive and useful. 
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